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The eigenfunctions of nested wells with incommensurate boundary geometry, in both hydrody- 
namic shallow water regime and quantum cases, are systematically and exhaustively studied in this 
letter. The boundary arrangement of the nested wells consist of polygonal ones, square or hexagonal, 
with a concentric immersed similar but rotated well or plateau. A rich taxonomy of wave patterns, 
such as quasicrystalline states, their crystalline rational approximants and some other exotic but 
well known tilings, is found in these mimicked experiments. To our best knowledge, these hydro- 
dynamic rational approximants are presented here for the first time in a hydrodynamic-quantum 
framework. The corresponding statistical nature of the energy level spacing distribution reflects this 
taxonomy by changing the spectral types. 

PACS numbers: 47.35.-i,47.54.-r,71.23.Ft, 73.21. Fg 



One would believe that a quasiperiodic wave pattern, 
that has an orientational order without periodic trans- 
lational symmetry 1], must be associated to an exter- 
nal single connected boundary. In such a case, classi- 
cal analogues which model features of quantum systems, 
and prove nontrivial properties of this last systems, have 
stirred interest. For example, acoustic Q and hydrody- 
namic quasicrystals Q were previously reported, and the 
existence of Bloch-like states has recently been proved in 
such systems 4]. In all of these results, there was an 
imposed global quasiperiodicity either by the boundary 
conditions or a dynamical source, which makes the ap- 
pearance of such patterns not so unexpected. However, it 
is also possible to confine quasicrystalline hydrodynamic 
modes within an inner isolated region of a bigger sur- 
face, as it was made in Ref. 5]. In that work, the ex- 
periment was realized under a linear regime and thus the 
fluid " sees" a double concentric non-connected boundary. 
Hence, it was essentially different to other reported non- 
linear quasicrystalline Faraday wave patterns 6] where 
the pattern is not conditioned by the shape of the bound- 
ary due to the non-linearity of the Faraday phenomenon. 

Although Ref. p provided new insights into the exper- 
imental generation of quasiperiodic patterns in the linear 
regimen, unfortunatelyit does not fulfills well the shallow 
water approximation Q, i-e., in which the surface wave- 



length is much larger than the liquid depth, making im- 
possible to state a suitable hydrodynamic-quantum anal- 
ogy. Since the quantum analogs of these confined hydro- 
dynamic modes may be of relevance to design quantum 
confinements echibiting quasiperiodic electronic states, 
or their rational approximants, here we systematically 
study the eigenfunctions corresponding to nested wells 
under the shallow water regime, and their corresponding 
quantum analogues. For this goal, we numerically solve 
the hydrodynamic equation, which in this case turns out 
to be the Helmholtz equation, and a thight binding ap- 
proach is applied for the quantum case. 

The nested wells of our study consist of a polygonal ex- 
ternal boundary, square or hexagonal, with a concentric 
similar but incommensurately rotated well or plateau. In 
the case of the squares, the rotation angle between both 
domains is 45°, and in the case of the hexagons, this an- 
gle is 30°. The bottom of the vessel was covered with 
a shallow liquid layer of depth hi, and the inner well or 
plateau with a depth /i2- 

In the hydrodynamic study we use the equation: 
d^rj{r,t) = V • [gh(r)Vrj(r,t)], Q where r] is the wave 
amplitude, g is the acceleration due to gravity and h(r) 
is the depth field, being the corresponding Helmholtz 



standing wave equations (V 



Vc?)¥(r)=0;i = l,2; 



where \I/(r) is the amplitude of the corresponding stand- 
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FIG. 1: Quantum octagonal quasicrystalline pattern (left) 
and its 3/2 rational approximant (right) in a double square 
nested well. Simulation parameters are as follows. For the 
quasicrystalline pattern, the well depth ratio is 1.16 and the 
level energy, respect to the 50Hz ground state, is —0.349174. 
For the rational approximant the well depth ratio is 1.1 and 
the level energy is —0.2246. 



ing wave, with the above mentioned adequate regions, 
1 (outer) and 2 (inner), and q are the phase velocities, 
uj/ki = (ghi) 1 / 2 , uj is the angular frequency and k is the 
wave number. 

We calculate the square nested well for a depth ratio 
h\/h,2 = 0.9,1.1,0.77, and 1.3 respectively, with about 
1700 eigenvalues and their corresponding eigenfunctions 
exhaustively checked for each case. The hexagonal nested 
well case is calculated for a depth ratio h\/h,2 = 0.9 and 
1.1, respectively, including about 1600 eigenvalues and 
their corresponding eigenfunctions in both spectra. Due 
to the high symmetry of the problem, a lot of doubly 
degenerate states appear. Any eigenfunction of spectra 
is available under request. 

Since in the present regime it is possible to establish 
an analogy with a quantum system, we also solve a tight- 
binding equation as described elsewhere [3, Q , with the 
same nested well geometry but providing it with a square 
atomic lattice. The correspondences with the hydrody- 
namic equations are the following: 

-Em,n — 4^f/l m?n / CI , 
tm,n,m — l,n = ^ra,n,ra+l,n = ^ra,n,m,n — 1 
= ^ra,n,m,n+l ^/l m?n / CI , 

where E m , n is the potential energy at the atomic site 
with integer lattice coordinates (m, n), /i m?n is the liquid 
depth at site (m,n), a 2 = L 2 /N 2 , where L is the length 
of the external polygon side and N is the discretization 
size (typically 201); and t m , n ,m',ri is the hopping integral 
between the nearest neighbour sites (m, n) and (m',n f ). 
Finally, in the present analogy, the dispersion relation 
corresponds to the equation: E = uo 2 a 2 / gh. Spectra with 
7.000 and 10.000 eigenvalues were respectively explored. 

As we shall see, we present here quasicrystalline wave 
patterns and their lower rational approximants linked to 
an hydrodynamic-quantum analogy problem. In Fig. ^ 
we show a quantum octagonal quasicrystalline pattern 
and its 3/2 rational approximant one. The first pat- 
tern is obtained under quantum boundary conditions, 
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FIG. 2: Hydrodynamic rational approximants of the octag- 
onal pattern: 3/2 (left top) and 7/5 (right top). The edge- 
lengths of the double square nested well are 80 and 35cm. In 
both cases the depth ratio is 0.9 and the frequencies for each 
case are 0.8376703Hz (3/2 approximant) and 0.9505245Hz 
(7/5 approximant). Below each pattern, their corresponding 
vibrational level distribution (VLD) are shown. 



i.e., Dirichlet conditions, and the second one, although 
it is also a quantum pattern, however it is obtained un- 
der hydrodynamic boundary conditions, i.e., Neumann 
conditions. Fig. ^ (left) conspicuously resembles to that 
pioneer experimental octagonal pattern early published 
in the very different non linear Faraday wave context 
by B. Christiansen, P. Alstrom and M. T. Levinsen 
(see Fig. 3 (bottom) of this reference). A similar qua- 
sicrystalline octagonal pattern as that shown in FigQ] 
(left) is obtained under our present hydrodynamic frame- 
work for a depth ratio of 0.9 and at an eigenfrequency of 
f = u/27r = 3.02632Hz. 

The quasicrystalline octagonal pattern is generated 
by using linear combinations of two square lattice 
vector bases, shifted between them by an angle of 
2 tan -1 (2 1 / 2 — l). By changing 2 1 / 2 for its respective 
rational approximant numbers obtained from the contin- 
uous fraction expansion 2 1 / 2 = [1;2], the approximant 
patterns are generated. These wave vectors were accu- 
rately obtained for each eigenfunction by means of its 
corresponding pattern Fourier analysis. 

In Fig|3 we also show both the hydrodynamic 3/2 
rational approximant pattern and the 7/5 one with 
their corresponding vibrational level distributions D(r]o) 
(VLD). The VLD is computationally obtained start- 
ing from the expression: D{t]q) = J _ q (dZ/|V77(r)|), 
where Q is the area of the rational approximant unit cell. 
By replacing the k-space by the real space, the above 
mentioned vibrational level distribution plays the same 
role in the corresponding eigenfunction as the well known 
density of states (DOS) in the whole spectrum [^.Hoj. 

In FigE3 we show an exotic pattern resembling the j3- 
Mn structure (or a-phase) which has been found coexist- 
ing as a crystalline approximant with the octagonal qua- 



3 




FIG. 3: Hydrodynamic pattern corresponding to the a phase 
obtained with the same geometrical parameters of the Fig- 
ure Eland with a depth ratio of 1.1 at a frequency equal to 
0. 7753178Hz. Its corresponding VLD is shown at the right. 



sicrystalline phase [ui] . As Fourier analysed in this work, 
this pattern grows spontaneously in this mimicked exper- 
iment with the same wave vectors of the 3/2 rational ap- 
proximant one but introducing the phase 2arctan(l/2) 
in a couple of non-orthogonal vectors of the correspond- 
ing 4D basis. The phasing out transforms the VLD of 
the pattern, as shown in Fig. |3f right), making it similar 
to that of the well known square lattice case 0, Hfl • It 
is remarkable that the present hydrodynamic experiment 
can give insight about a basic crystallographic structural 
problem. To the authors best knowledge, this is the first 
time that the a-phase is described in the density wave 
framework. 

In a similar scenario, Fig^I shows two patterns gen- 
erated in a hexagonal nested well and their correspond- 
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FIG. 4: Hydrodynamic triangular pattern (left) and the do- 
decagonal 19/11 rational approximant (right). The edge- 
lengths of the double hexagonal nested well are 40 and 17.5cm 
respectively. The depth ratio and frequency for each case are: 
1.1, 0.8155099Hz (triangular) and 0.9, 2.2416974Hz (19/11 
approximant). Below each pattern, their corresponding VLD 
are shown. 
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FIG. 5: Poisson distribution of nearest level spacing in the 
whole level spectrum corresponding to the square symmetric 
nested wells case. Bars are experimental results and dashed 
line is the theoretical distribution. 



ing level distributions. The quasicrystalline dodecagonal 
wave pattern can be also generated using two ternary 
wave vector sets shifted between them by an angle of 
2 tan -1 (2 — 3 1 / 2 ). Their rational approximant patterns 
appear when 3 1 / 2 transforms to its corresponding ratio- 
nal approximant numbers. These rational numbers are 
obtained starting from the continuous fraction expansion 
3 1 / 2 = [1;T72]. FigH shows the 2 and 19/11 above men- 
tioned rational approximant wave patterns. The VLD 
of the first pattern, or triangular pattern, is well known 
0, ^J. The spectral problem and the rich taxonomy 
of eigenfunctions of this double-hexagonal experiment is 
similar to that previously described for the double square 
well. 

As is well known |l2|, there is a relationship between 
the geometry of the boundary, the eigenfunction nature 
and the statistical properties of the spectrum. An in- 
tegrable system, such as a polygonal billiard pjj, has 
a statistical distribution of energy spectral fluctuations 
(P(«s)), where s is the distance between nearest energy 
levels, that fits the Poisson distribution. In Fig|SJ we 
show that this is case for the P(s) of the system under the 
symmetric configuration of the nested wells. If the sym- 
metry of such system is broken, by slightly rotating and 
shifting the inner well, then all the inherent degenera- 
cies are removed as it was made in a similar early study 
14]. In such a case, the eigenfunctions show "scarred" 
localized surface waves [l5j as shown in Fig|Sl and P(s) 
becomes a semi-Poisson distribution. 

We have shown that a geometrical arrangement of two 
nested wells, each of them with polygonal symmetry, can 
give a very rich and complex behavior even if a sim- 
ple linear differential equation, valid for hydrodynamics 
and quantum mechanics, is used as physical description. 
Quasiperiodic, approximant and periodic wave patterns 
have been obtained in this work. The present approach 
gave a new insight on how to generate certain rational 
approximant s in the density wave framework. The ob- 
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served spectral fluctuations are consistent with the ob- 
tained wave patterns in both symmetrical and chaotic 
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FIG. 6: Scarred chaotic wave pattern (top) obtained when the 
degeneracies are removed by rotating and shifting the inner 
well. An example is shown corresponding to a depth ratio 
of 0.77 and at a eigenfrequency of 4.14917Hz (bottom), and 
its Semi-Poisson distribution of energy spectral fluctuations. 
Continuous line represents the theoretical distribution and 
bars are experimental results. 
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